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Abstract 

The possible phases and the nano-scale particle correlations of two-dimensional interacting dipo¬ 
lar particles is a long-sought problem in many-body physics. Here we observe a spontaneous 
condensation of trapped two-dimensional dipolar excitons with internal spin degrees of freedom 
from an interacting gas into a high density, closely packed liquid state made mostly of dark dipoles. 
Another phase transition, into a bright, highly repulsive plasma is observed at even higher exci¬ 
tation powers. The dark liquid state is formed below a critical temperature D ~ 4.8iC, and it is 
manifested by a clear spontaneous spatial condensation to a smaller and denser cloud, suggesting 
an attractive part to the interaction which goes beyond the purely repulsive dipole-dipole forces. 
Contributions from quantum mechanical fluctuations are expected to be significant in this strongly 
correlated, long living dark liquid. This is a new example of a two-dimensional atomic-like interact¬ 
ing dipolar quantum liquid, but where the coupling of light to its internal spin degrees of freedom 
plays a crucial role in the dynamical formation and the nature of resulting ground state. 
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Cold two-dimensional dipolar gases are an intriguing and theoretically challenging ex¬ 
ample of many-body interacting quantum fluids. Theoretical attempts to understand the 
phases and correlations of such dipolar fluids revealed a very rich behavior, with predictions 
of several exotic quantum and classical phases and correlation regimes [1-7] that are absent 
in weakly interacting Bose gases. These predictions sparked intensive experimental efforts 
in both atomic and molecular systems [8-12], as well as in two dimensional (2D) dipolar 
excitons in semiconductor bilayer systems (also known as indirect excitons - IXs) [13-16]. 
These IXs are composed of a coulomb-bounded, yet spatially separated electron-hole pairs in 
an electrically biased semiconductor double quantum well (DQW) structure [17, 18]. Recent 
theoretical analysis of the correlations in such a dipolar IX system predicts that in a very 
large range of densities the onset of quantum degeneracy is accompanied by the set-in of 
multi-particle interactions and correlations, so that dipolar IXs should behave as a strongly- 
interacting liquid rather than as a weakly interacting Bose-Einstein condensate (BEG) [5, 6]. 
What makes this IX system even more unique is the additional internal spin degrees of free¬ 
dom of excitons in GaAs-based DQWs, allowing to explore the effect of spin on quantum 
correlations. The IX can be found in four different spin states, namely S = ±1,±2. The 
IXs with S' = ±1 are optically active and are named ’’bright” while the ones with S = ±2 
are optically inactive and are named ’’dark” [19]. This dark IX thus have a much longer 
lifetime. Due to spin-dependent exchange interaction between electrons and holes, the dark 
IXs are slightly lower in energy than the bright IXs [20], and therefore it was predicted that 
the ground state of the dipolar exciton fluid is dark and that a macroscopic phase transition 
to a dark BEG should be observed [21]. 

Erom the theoretical point of view, it seems that the picture is still incomplete, as no 
theory takes into account both dipolar interactions and the internal spin degrees of freedom 
to predict the origin and nature of the ground state of the IX fluid. On the experimental side, 
there are reports on the onset of coherent phenomena expected from a quantum degenerate 
dilute gas of IXs [14, 22-24], but other works show experimental evidence for a formation 
of a highly correlated interacting fluid [15] and for an exciton liquid [25] and of a notable 
spontaneous darkening of a correlated dipolar exciton fluid below a certain temperature 
[15]. It is therefore an open challenge to fully and unambiguously determine the formation 
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mechanism and properties of the low temperature highly correlated collective dipolar state. 
Can the combination of the strong dipolar interactions, quantum degeneracy, together with 
the very long lifetime of the dark particles lead to a long lived, dark quantum dipolar liquidl 

In this paper we report on the first observation of a spontaneous gas-liquid-plasma tran¬ 
sitions of a trapped IX cloud at low temperatures, and show that the liquid phase is made 
up of mostly condensed dark particles, that it has the characteristic parameter range of 
an incompressible closely-packed quantum liquid, and that the underlying physics probably 
does not result from purely dipolar forces, but from more complex interactions at high den¬ 
sities, which have to do with the short range interactions due to the internal spin degrees of 
freedom. 

To have a well defined and highly controlled system, we trap a cloud of IXs, optically 
excited in a GaAs based DQW structure under an electrostatic trap gate [26], as is depicted 
schematically in Fig.l. The IXs are excited non-resonantly at the center of this quasi- 
parabolic electrostatic trap (denoted Xtrap), and its photoluminescence (PL) is measured 
under steady-state conditions (see Fig.2(a,b) for the Xtrap geometry and its energy profile 
respectively). 

The spatial-spectral emission profile of the IXs is imaged on a CCD, allowing us to follow 
the size and shape of the IX cloud. Furthermore, we use a constant energy line method 
(CEL, see [27] and SI) in which we can keep the energetic and electric environment of the 
exciton fluid fixed. This is achieved by adjusting the magnitude of the applied electric field 
F on the DQW to keep the band tilting and thus the difference between the direct exciton 
(DX) and the IX transition energies, Edx — Ejx fixed while varying the optical excitation 
power and temperature. With this method we can extract information on the relations 
between the energy resulting from the many-body interactions between the IXs - AE [15] 
and their density, while fixing the single-particle properties of an IX (either bright or dark). 
This is because these single IX properties, such as its radiative and non-radiative lifetimes 
and its effective dipole size d are determined by the same band tilting which also uniquely 
determines Edx — Ejx- 

Four exemplary spatially resolved PL (along a central cross-section of the circular Xtrap) 
are shown in Fig.2(c-f) for non-resonant excitation at the center of the Xtrap. These measure- 
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(a) 





FIG. 1. (a) A sketch of the sample structure, showing two quantum wells (light blue) separated 
by a thin barrier (red). An electric field is applied perpendicular to the growth axis between a 
transparent lOnmTi top circular gate and doped GaAs layer at the bottom. A non-resonant 
optical excitation using a GW diode laser creates a trapped IX fluid under the gate, as illustrated. 
The tunneling current through the sample is also monitored, (b) Schematics of a calculated energy 
bands diagram for the structure, showing the conduction and valence bands and the electron and 
hole wave functions, respectively. IX marks the indirect exciton transition, whereas DX marks 
the direct exciton. (c) The first energy levels of direct exciton transitions (dashed lines), of either 
electrons (e) and heavy holes {hh) or light holes {Ih) are marked, as well as the laser excitation at a 
wavelength of 780nm (solid arrow), illustrating the non resonant excitation of the quantum wells. 

ments show two different temperatures and excitation powers, while the applied voltage was 
adjusted in each of the measurements to maintain a constant energy {E^x — Ejx ~ 13meV) 
for all of them. Two excitonic PL lines are observed for all the spectra, the high energy 
line (very weak) corresponds to the DX while the low energy line correspond to the IX. The 
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FIG. 2. (a) A microscope image of a 20/im diameter electrostatic trap (Xtrap) surrounded by a 
guard gate. The bright spot at the center is the exciting laser image, (b) Spatial-spectral profile 
of the Xtrap, obtained by scanning a weak and focused excitation point along the line marked 
in (a), (c-f) Spatial-spectral PL intensity for temperatures T = 1.9,6.2iF and excitation powers 
P = 29,1891nVF, all taken on a CEL by adjusting the voltage. The horizontal axis shows the PL 
energy, and the vertical axis is a cross-section of the position along the Xtrap. The peak with the 
lower energy {Ejx ~ 1.541eE) is PL from the indirect exciton (IX), and the weak emission with 
the higher energy {E^x ^ 1.554eE) is PL from the direct exciton (DX). Right: spatial intensity 
profiles of (c-f), obtained by integration along the spectral axis, showing (c,d) very minor spatial 
expansion of the PL for increasing P at low T, in contrast to (e,f) a substantial increase of the 
PL area with excitation power at high T. Bottom: the extracted spectra from (c-f), obtained by 
integration along the spatial axis, illustrating (c,e) a narrower line-shape at low T and, (d,f) a 
reduction of the linewidth to a symmetric and narrow emission line-shape as the excitation power 
increases. 
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PL is then spectrally integrated to obtain spatial profiles of the emission (Fig.2 right) and 
spatially integrated to obtain a single spectrum for each measurement (Fig.2 bottom). Such 
spectra are obtained and analyzed for different experimental conditions as is described next. 

Figure 3(a) shows the total intensity Iix{T) extracted from the area under the spectral 
line of the IXs, obtained at different temperatures for a fixed excitation power (P 308nbF) 
on a CEL with a fixed Ejx ~ 1.541eE. Also shown is the applied electric field F that 
was required in order to maintain the CEL at the aforementioned energy. By performing 
such CEL experiments for excitation powers P spanning almost two orders of magnitude 
for each T, we extract the many-body interaction term AE{T, P) in a straightforward way: 
AE{T,P) = eF{T,P)d — eF{T,0)d where F{T,P) is the applied field required for a fixed 
Edx — Ejx at a given (T, P), and d is the effective distance between the centers of the 
two QWs (see SI for more detailed information on this method). These AE values and the 
spectral linewidth of the IXs PL are plotted as a function of T for different P values in 
Fig.3(b,c). 

In the higher temperature range (T > 4.8iF), AE shown in Fig.3(b) decreases with 
decreasing temperature but increases with increasing excitation power. The increase of AE 
with P is not surprising: the magnitude of the repulsive dipole-dipole interactions increase 
with increasing total density n, i.e., AE — C{n,T) ■ n (where C is the interaction-induced 
particle correlation function [5]), and n increases with increasing P. The decrease of AE 
with decreasing T (for a fixed P) at this temperature range was also predicted theoretically 
[5, 6, 28] and observed experimentally [15], and it arises from increased interaction-induced 
particle correlations as the IXs become colder, resulting with a reduction of C and thus 
of AE . Already at this higher temperature range a reduction of the IX linewidth with 
increasing P for all temperatures above 4.8^ is observed, as seen in Fig.3(c). Such a strong 
linewidth reduction with increasing particle density indicates that the increase of density 
reduces the fluctuations of the interaction energy around its average value. This density- 
induced linewidth narrowing is a strong indication for multi-particle interactions and for 
deviation from a gas behavior [25]. 

Surprisingly, moving to the low-T range (T < 4.8P), AE changes its trend and sharply 
rises with decreasing T, as can be seen in Fig.3(b). AE slightly ’’overshoots” to a maximum 
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FIG. 3. Results obtained for a CEL with Ejx = l.SdleE. (a) Total intensity of the IX line - 
Iix (left axis, blue circles) and the electric field F (right axis, orange triangles) needed in order to 
maintain this CEL, as a function of temperature for a fixed excitation power of P ^ 308nVF. (b) 
The blue-shift energy AE and (c) the IX PL linewidth as a function of temperature T for several 
different excitation powers. In the high temperature regime (T > Tc), AE decreases with tem¬ 
perature, indicating increased particle correlations. For T < Tc^ AE increases as the temperature 
is lowered and saturates to an almost fixed value, regardless of the excitation power. This marks 
the condensation of IXs into a fixed, high density liquid state. The reduction of the PL linewidth 
with decreasing T to a small fixed value points at reduced density fluctuations, as expected from 
a liquid phase. The inset in (c) shows the dependence of the PL linewidth on excitation power for 
three different temperatures. 
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at T ~ 2K and then slightly decreases and reaches a constant value, AEi ~ 10.5mel/, at the 
lowest measured temperature, for all excitation powers P. The PL linewidth also converge 
below T = 4.8iL to a constant low value at the lowest measured temperature, as can be seen 
in Fig.3(c). At this temperature, the linewidth also becomes independent of the excitation 
power, as is seen in the inset. The sharp increase of AE below 4.8iF indicates an increase of 
the total dipolar fluid density n [29]. Furthermore, the fact that AE —> AEi at the lowest 
measured T, together with a reduction of the PL linewidth to a small, fixed value suggests 
that the trapped fluid reaches a fixed density ni regardless of the excitation power. This 
observation is consistent with a transition into a condensed liquid state, that has a high 
density, fixed by the constraints of the physical system. 

Perhaps the most striking evidence for the condensed phase formation and its properties 
can be seen in Fig.4. Since in Fig.3 the condensation at low temperatures was observed for 
all the measured excitation powers, we expanded the excitation power range even further 
to cover five orders of magnitude, and capture how the phase transitions are driven by the 
increased rate of optical excitation. The top-right panels of Fig.4 show the spatial-spectral 
emission profile of the IXs at T = l.bK < Tc- Initially, increasing the power results in an 
expansion of the IX cloud in the trap, as is expected from the mutual repulsive interactions 
between IXs. However, above a certain threshold power, the cloud shrinks with increasing 
excitation power and reaches a fixed size which almost does not vary over about two orders 
of magnitude of excitation powers. This marks the condensation to the fixed density liquid 
phase. When the power is further increased above another threshold value, the cloud starts 
to expand again in a fast rate. This fast expansion is a signature of a transition to an e- 
h plasma, driven outwards by the intra-layer e-e and h-h Coulomb repulsion. The spatial 
profile FWHM and the corresponding AE dependence on the excitation power at T = l.bK 
are shown by the blue symbols in Fig.4(a,b) respectively. It can be seen that the initial low 
power expansion turns into a very strong contraction that is accompanied by a large increase 
of AE, indicating a spontaneous large increase of density (which we will discuss later on). 
This is followed by an intermediate power regime (corresponding exactly to the two orders 
of magnitude power range of Fig. 3), in which both the spatial FWHM and AE are nearly 
constant, corresponding to an almost fixed density, regardless of excitation power. The high 


power regime shows a sudden and sharp increase of both the spatial FWHM and of AE, 
consistent with the picture of a strongly driven indirect e-h plasma. The orange symbols 
show the behavior at T = bAK > Tg, above the condensation temperature. Indeed, non of 
the features of the condensation are observed, and the IXs expand and increase their density 
monotonically, consistent with a repulsive dipolar gas phase. 

We can now estimate the saturated density of the condensed phase, ni, from the saturated 
value AEi. Such an estimate requires a knowledge of the spatial particle correlations and 
C. In Ref. [6] a theoretical estimate for C was given for a dipolar exciton liquid. The 
density estimated from this liquid model yields ni ~ (by setting 
the experimental value AEi = lO.bmeV and the dielectric constant k = 13. A more self- 
consistent estimate, taking into account the kinetic energy from quantum motion (see SI), 
gives u; ~ 9 • This density corresponds to an average inter-particle distance 

(t ^ ^ Ttj^ ^ 30nm. This average inter particle distance is smaller than the effective 

/ 2 j2 \ 

interaction radius of the dipolar exciton [6] ro = ( j ~ bbnm at T = 2K, which 
means that at this density the fluid is deeply inside the multi-particle interaction regime, 
thus the picture of a dipolar liquid state is self-consistent. 

Remarkably, (r) is between twice to three times the exciton Bohr radius (ax ~ 10—15nm), 
namely {r)/ax ~ 2 — 3. Therefore, at this density the IXs are almost closely packed and 
any further increase in the density significantly increases the overlap of the excitons wave 
functions. Such an increased overlap should result in a large deviations from the pure dipolar 
repulsion due to short-range direct and exchange coulomb terms. This specific inter-particle 
distance might point to the physical mechanism which both drives the spontaneous spatial 
shrinkage and determines the almost fixed density and thus the low compressibility of the 
condensed liquid. Recent numerical simulations of the two-body interactions of dipolar 
excitons have suggested that the interaction of two dipolar excitons can become slightly 
attractive as the inter-particle distances approach (r) ~ 2 — 3ax , and then strongly repulsive 
scaling as 1/r® at even shorter distances [30]. This inter-particle distance is similar to the 
estimate of the current experiment. 

A natural question now is what can be the mechanism responsible for the observed sharp 
increase in the IX cloud density as the excitation power increases above the critical power? 
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FIG. 4. Results obtained for a CEL with Ejx — 1.538eE with excitation power ranging almost 
5 orders of magnitude. Right panels: spatial-spectral photoluminescence images of the trapped 
IX fluid taken at T = l.bK with increasing excitation power. Left panels: (a) AE and (b) 
spatial FWHM of the exciton fluid as a function of excitation power at T = l.bK (blue curves), 
compared with those measured at T = bAK (red curves). A transition between three different 
regimes as the excitation power increases is observed only for the low-temperature experiment, (c) 
Estimated density of the exciton cloud ritotal a function of excitation power at T = l.bK. The 
data obtained from photon-counting (circles) is presented as twice the bright exciton density, 2ni, 
(assuming = n^), together with an expected density range (triangles, shaded region) that is 
calculated from the measured blueshift energy AE (see text for details). The black vertical arrow 
indicates the difference ritotal ~ ‘^'^h — ~ (d) Estimated values of length scales for an exciton 

liquid: ax ~ 15nm is the exciton Bohr radius, tq ^ lOnm is the quantum mechanical wavefunction 
radius of an exciton, tq ^ bbnm is the classical minimal distance between excitons and At ~ 200nm 
is the thermal de-Broglie wavelength (both calculated at T = 2K). See text and SI for details on 
the derivation of those values. 
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In principle, a sharp increase in the total density under a fixed excitation power can result 
from a sharp increase in the radiative lifetime of the particles. Noting however that the whole 
experiment was done under single CEL, where the radiative lifetime of a single IX is kept 
constant, such an increase is not plausible. The other option is thus that the thermal equality 
between the populations of dark and bright IXs breaks, and that there is a spontaneous 
macroscopic accumulation of dark IXs above the critical excitation power. As dark IXs are 
optically inactive and have much longer lifetimes [15], this should result in a sharp increase 
of the total density for a fixed excitation power, as is indeed observed in the experiment. 

To check this latter picture, we extracted the total particle density, n = 71^, + rid, from the 
interaction energy AE, as plotted in Fig.4(c). To avoid any model assumption about the 
relation between n and AE through the unknown correlation parameter C, the shaded blue 
region covers all possible values of C, where the lower boundary is given by the uncorrelated 
gas limit (mean-field) and the upper boundary by the highest correlated limit of a dipolar 
crystal [5, 18]. The sharp increase of n at the critical power of the gas-liquid phase transition 
and its saturation up to the second critical power are clear. This behavior of n extracted 
from AE is then compared with the total density estimated from the bright IX density ni, 
alone, which we very carefully extracted from the intensity of IX emission (see SI for details). 
If we assume a thermal equality between bright and dark IX populations, i.e., rid = then 
n = 2nft. This total density is plotted by the blue circles in Fig. 4(c). As can be seen, 
the density of bright IX smoothly increases with excitation power and does not exhibit any 
sharp increase at the phase transition, and thus it gives a much lower density than the density 
inferred from AE. Since AE inherently depends on the total density n, while an assumption 
of thermal equality of bright and dark IXs population was used to extract n from the PL 
intensity, it is clear that this assumption strongly underestimates the total density all the 
way to the second critical power where the transition to a plasma is observed. As such a 
strong underestimate is not observed for the data above Tc (see SI), this strongly supports 
the picture of condensation of the IXs population into a dark state [21]. It is interesting to 
note that the thermal de-Broglie wavelength of excitons at this temperature is estimated to 
be At = h/\/2Mxk^ ~ 200nm. Therefore, (r) -C At so that the IX cloud is quantum 
degenerate and this should support the condensation of dark IXs [31]. 
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Finally, another important length scale that characterizes the liquid phase the is the 
quantum radius of the IX. This quantum radius, tq, is the ground-state wavefunction radius 
of an exciton confined in the parabolic-like potential resulting from the interaction with 
its surrounding excitons in the liquid. With a density ni given above, we can estimate 
vq ~ lOnm (see SI and Ref. [6] for the calculation of tq). Such a comparison yields {r)/2rQ ~ 
1.5. This result suggests that the quantum mechanical motion of the dipoles due to the 
uncertainty principle (their zero-point motion) is of the order of the inter-particle distance. 
Two important properties of the liquid can be inferred from this ratio: the first is that 
quantum mechanical fluctuations prevent formation of long range order and crystallization. 
The second is that it might suggest that the IX liquid has strong quantum mechanical 
properties, similarly to liquid for example. Indeed, the same calculation (see SI) shows 
that the kinetic energy due to the quantum motion Eq is much larger than the classical 
kinetic energy ksT, and that it contributes about 25% of the total liquid energy. 

In conclusion, we observed phase transitions of a dipolar exciton fluid from a gas to a 
highly correlated dark liquid, and then to an e-h plasma, as the excitation power is increased 
at low enough temperatures. The liquid state, which results from a condensation of dark IX 
particles, has a well defined density, indicating a significant reduction in its compressibility. 
The analysis of the liquid’s energy and other properties suggests that the inter-particle 
distance approximately corresponds to a close-packing of the dipolar excitons, and that the 
liquid might have a significant quantum uncertainty energy contribution, quite similar to 
liquid '^He. This is a new and exciting example of a collective state of a two-dimensional 
interacting system with internal spin degrees of freedom. Furthermore, this artificial atomic- 
like fluid resides inside a low-dimensional semiconductor structure, and it is a unique example 
where coupling of light to internal spin transitions plays a crucial role in the formation 
dynamics and in determining the collective ground state properties. 
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Supplementary Information 


SUPPLEMENTARY NOTE 1 - RAW DATA FOR DIFFERENT TEMPERA¬ 
TURES AND EXCITATION POWERS 

Figure SI shows some raw data of photoluminescence, obtained at different temperatures 
and excitation powers, all for a constant energy line Eix ~ 1.541eU. One can see a very 
minor increase along the Xtrap spatial cross-section for the highest excitation power at 
T = 1.5K, which translates into a slightly increased trap area by a factor of less than 1.4. 


T=1.5K, P= 29nW T=1.5K, P= 308nW T=1.5K, P=1891nW Spatial profile, T=1.5K 





T=4.8K, P= 29nW T=4.8K, P= 308nW T=4.8K, P=1891nW Spatial profile, T=4.8K 





T=6.2K, P= 29nW T=6.2K, P= 308nW T=6.2K, P=1891nW 


Spatial profile, T=6.2K 
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FIG. SI. Spatial-spectral photoluminescence images recorded for temperatures T ^ 1.5,4.8, 6.2K 
and excitation powers P 29, 308,1891nVF, all for constant energy line Ejx ~ 1.541ey. The 
column to the right shows the corresponding spatial profiles for each temperature, obtained by 
integration along the spectral axis, and the legend shows FWHM of the curves (color ordering is 
for low, medium and high power). 
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SUPPLEMENTARY NOTE 2 - EXTRACTION OF THE BLUE-SHIFT ENERGY 


AE FROM THE APPLIED FIELD 


The electric field F that is needed in order to maintain a specific constant energy line 
(i.e. a fixed Ejx) is calculated using F = V/L, where V is the applied voltage and L is the 
active sample thickness. Such data is shown in Fig. S2 as a function of temperatures for 
several different excitation powers. 
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FIG. S2. The electric field F that is needed in order to maintain a specific constant energy line of 
Eix ~ 1.541eF as a function of temperature, for several different excitation powers. 

In order to extract the blue shift energy AE(T, P) from F(T, P), we note that 

Eix = Edx - eF(T, P)d + AE{T, P). (SI) 

where the direct exciton energy Edx has a very weak dependent on the density, temperature 
and the applied field (quadratic stark shift) in the measured parameters range. Now, in the 
limit of a very weak excitation power the indirect exciton gas in the trap becomes very dilute, 
and the blue shift energy is negligible: AE —?• 0 so one gets 

Eix = Edx — eF{T,0)d (S2) 

For a fixed Edx — Eix (he. on a given CEL) we can combine the above equations and get 
that 

AE{T, P) = eF{T, P)d - eF{T, 0)d. (S3) 
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Extrapolation of F vs. P @ T=6.2 



FIG. S3. The electric field F as a function of excitation power P, recorded at a temperature 
T m 6.2K. Blue circles indicate measured data points while the red curve is an interpolation 
(between points) or extrapolation (towards P — 0) of the data. 

Next, we extrapolate the data of electric field F vs. excitation power P towards zero to 
extract F{T, 0) in the limit of zero excitation power, as is shown in Fig. S3. 

The temperature dependence of F(T, 0) is expected to be negligible since the band gap 
dependence with T is essentially flat in this temperature range and there are no multi¬ 
particle interactions in this limiting low density. Therefore we chose the data of T ~ Q.2K 
for the extrapolation process and the same electric held F{T = 6.2K,0) is used to deduce 
the blue shift energy for all temperatures. This temperature is chosen since for (T > Tc) 
there is a clear power dependence of F, and no accumulation of dark particles is expected. 
The extrapolation of F{T,P) towards F(T, 0) for other temperatures above T > 4.8iF gives 
similar results within less than 5% error. 

SUPPLEMENTARY NOTE 3 - INTENSITY AND BLUE-SHIFT ENERGY FOR 
DIFFERENT CONSTANT ENERGY LINES (CELS) 

The results presented in the paper are not restricted to a specific choice of indirect exci- 
ton energy Fix, but rather show similar features for the whole range of energies that was 
measured. Figure S4(a) shows the recorded PL intensity as a function of temperature for 
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different Ejx values (at fixed excitation power). Figure S4(b) shows the applied electric field 
F needed to maintain the CEL as a function of temperature. As one can expect, lower Ejx 
requires stronger electric field in agreement with the linear Stark effect. We can then extract 



123456789 10 

T(K) 

FIG. S4. (a) Measured PL intensity as a function of temperature for different CEL. All curves 
are for a fixed excitation power P — 308nIF. (b) The applied electric field F needed in order to 
maintain the CEL as a function of temperature, for the same set of measurements as in (a). 

the blue-shift energy AE as is described in Supplementary Note 2. These results are shown 
in Figure S5, and display essentially the same features as shown in the main text. 


SUPPLEMENTARY NOTE 4 - PHOTOLUMINESCENCE INTENSITY VS. TUN¬ 
NELING CURRENT 

While maintaining a constant energy line of Ejx = 1.541eU, we recorded the total IX 
photoluminescence (PL) and tunneling current as a function of temperature for different 
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FIG. S5. The blue shift energy AE as a function of temperature, for several different excitation 
powers. Results are shown for two CELs with Ejx = 1.535eE and Ejx = 1.538eE. 

excitation powers. The results are shown in Fig. S6, where we plot the PL intensity and 
tunneling current normalized by the excitation power. One can see an increasing PL intensity 
as the temperature is lowered down to about T 3K, where the intensity begins to decrease 
(darken). The tunneling current shows roughly an opposite pattern (except for the two lowest 
excitation powers, where the actual values are within or below our measurement resolution). 
These findings can be explained in fight of the competing nature of radiative and non- 
radiative decay mechanisms. At high temperatures, there is a lower fraction of excitons 
occupying the light cone, and thus non-radiative mechanisms such as tunneling current are 
more dominant. As we cool down, more excitons enter the light cone and become radiative, 
so the intensity increases and the current decreases. 

However, there is a clear change of trend below T 3K, from increasing to decreasing 
intensity as we lower the temperature. This suggests that radiative processes become less 
probable, as more excitons reside in dark states, and are thus compensated by an increasing 
tunneling current. To check this hypothesis, we carried out a careful and detailed analysis of 
the data which allowed us to compare the total number of emitted photons to that of electrons 
tunneling through the sample (for details see supplementary note 5 next). The results are 
presented in Fig. S7 and indeed show a very good agreement between the decrease of the 
photons’ rate at the lower temperature range and the increase of charge carriers rate. This 
indicates that the ” loss” of photons is accounted for by the increase of e-h pairs in the current 
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FIG. S6. (a) The total IX intensity Iix and (b) the DC current flowing through the sample, as 
a function of temperature for several different excitation powers. The curves in are normalized to 
their corresponding laser power P. Error bars are shown only for a single power level. 

(each exciton contributes one e-h pair if it tunnels out non-radiatively), so conservation of 
particles is maintained. 


SUPPLEMENTARY NOTE 5 - BRIGHT EXCITONS DENSITY ESTIMATION 


In order to calculate the bright exciton density Ub from the measured photoluminescence, 
we assume a thermal distribution of excitons with uniform density across the trap area. In 
this case, we can write a rate equation for Ub of the following form: 

dnb /3nb 


— cTif^ T G 


(S4) 


dt T 

where (3 = Urad/nb is the fraction of bright excitons occupying the light cone (/5 is 
temperature-dependent) [15], r is the single exciton radiative lifetime, G is exciton gen- 
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FIG. S7. (a) The total emitted photons rate and (b) the charge carriers flowing through the sample 
rate , as a function of temperature for several different excitation powers. Color code is the same 
as in Fig. S6. 


eration rate and c is a bi-molecular recombination coefficient describing radiative emission 
due to two-particles processes (where cn^ <C /3/r). In addition, due to the limited numerical 
aperture of our experimental apparatus, only a relatively small fraction of photons that are 
emitted from recombining bright excitons are collected and reach the detector (CCD). We 
denote this fraction by 7 == ricoiiectedlnb. 

In general, these ratios depend on the distribution of excitons as a function of energy 
since they arise from conservation of energy and in-plane momentum (see Fig. S 8 ). Figure 
S9 shows for example 7 for the case of a uniform emission distribution (which means no in¬ 
plane momentum conservation) compared with a thermal Maxwell-Boltzmann distribution 
of excitons with momentum conservation. The error bars in the following estimations and 
those presented in Fig. 4(c) of the main text account for both limiting cases. 
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FIG. S 8 . Left: an illustration of a recombining exciton. k\\ is the IX’s in-plane momentum, while 
and k are the photon’s perpendicular and total momentum, respectively. Right: schematics of 
a thermal distribution of bright IXs kinetic energy. Excitons with energy below some value E* are 
inside the light cone and thus radiative, while photons originating from excitons with energy below 
are collected into the numerical aperture of our system. 



FIG. S9. The fraction of excitons whose emitted light is collected by the experimental system 
numerical aperture of NA=0.42, 7 = assuming a Maxwell-Boltzmann thermal distribution 

of excitons with momentum conservation (solid, blue), compared with a uniform emission pattern 
into a 47 r solid angle (dashed, red). 

The intensity / (counts/unit time) recorded by our detector is given by 

I = r) ■ A (^^ + cn^) (S5) 
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where A is the exciton cloud area which is imaged onto the CCD, and rj is the combined 
loss of our experimental setup and the conversion efficiency of photons into counts. After 
measuring and accounting for all losses incurred by the optical elements of our setup and 
carefully calibrating the photons/count ratio in our system, we can extract the bright exciton 
density, r was measured independently using a time resolved measurement described in detail 
in Ref. [15] so the only fitting parameter (i.e. the unknown parameter) in our model is c, 
which is determined by the best agreement between the total density as extracted from AB 
and that of 2ni, for T = 5AK, as shown in Fig. SIO. We note however that the fit weakly 
depends on c, since assuming c = 0 only slightly changes the quality of the fit. 



10 "^ 10 ‘‘ 
Excitation power (W) 


FIG. SIO. Estimated density of the exciton cloud njota/ as a function of excitation power at T = 
5AK and Ejx — 1.538eF. Twice the bright exciton density, 2nft, obtained from photon-counting 
(circles) is presented together with the density range (triangles, shaded region) that is calculated 
from the measured blueshift energy AE limited by an uncorrelated mean-field interactions (lower 
limit) and the highest correlated state of an excitonic crystal (upper limit). 


SUPPLEMENTARY NOTE 6 - DENSITY ESTIMATION USING AN INTER¬ 
ACTING QUANTUM LIQUID EXPRESSION 

In the main text, we used the theoretical classical liquid interaction energy adopted from 
Ref. [6] (see Eq. 5.5 of this reference) to deduce the exciton density. However, this expression 


23 








neglects the contribution from quantum confinement energy that comes from ’’caging” the 
particle and thus adding a zero-point energy term which becomes not negligible in this 
density range. In order to have a more accurate density estimation, we can add the quantum 
correction to the interaction energy (Eq. 5.2 of Ref. [6]) to have: 


AE 


lOe^d? 


K 






kM 


-n 


5/4 




Ec 


Eq 


(S6) 


Solving this equation numerically for a value of AE = 10.5meE, one gets a density n ~ 
9 • Plugging back this density into the equation yields Ec = 7.7meV and Eq = 

2.8meV, so indeed the kinetic quantum confinement energy is significant. 


SUPPLEMENTARY NOTE 7 - SPATIAL COHERENCE MEASUREMENTS 

We performed spatial coherence measurements on the same IX system discussed in this 
work, covering a wide range of excitation powers and temperatures, using a flip-image in¬ 
terferometry [32], and were unable to detect any spatial coherence of the PL over distances 
larger than the spatial resolution of the measurement system Ifim), even at excitation 
powers and temperatures where the anomalous behavior reported in this work was found. 
We used a Michelson flip-image interferometer, shown in Fig. Sll, to measure the spatial 
coherence of PL emitted from indirect excitons in the trap. The spatial resolution of the 
optical system as given by A/2NA 1pm for our parameters A ~ SOOnm and NA = 0.42. 

The experimental method is as follows: 

• Aligning the two arms of the interferometer to a sub-picosecond time delay using an 
ultra fast laser source (to ensure coherence time longer than the reciprocal of the PL 
linewidth). 

• Passing the PL from the sample through this configuration. 

• Scanning the phase of one arm for few times 27r while recording images, and detecting 
the oscillation frequency in each pixel due to the movement of the fringes. 
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FIG. Sll. A Michelson flip-image interferometer. 

• Using Fourier analysis to obtain the amplitude of oscillation in the detected frequency 
which, after background substruction and proper normalization, represents the degree 
of coherence. 

The above procedure was performed over a range of experimental conditions similar to 
those shown in the main text. Nevertheless, no extended spatial coherence was detected, as 
illustrated in Fig. S12. This lack of coherence can be consistent with a formation of a dark 
liquid, where even if the coherence is not strongly suppressed by the strong interactions, it is 
mostly manifested in the dark exciton part and therefore no spatial coherence is observed in 
the bright, emitting part. We note that other experiment on cold IX fluids were also unable 
to detect any spatial coherence down to very low fluid temperatures [25, 33]. 
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FIG. S12. A spatial degree of coherence map of the PL from trapped IX. Experimental conditions 
are T = 1.7K and P = 2/iIF. The dashed line marks the boundaries of the Xtrap. 
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